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Abstract 
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duction and decay of the delta resonance. One of our results is the agreement with the MiniBooNE 
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muon and electron type neutrinos, where effects of the lepton masses are visible. 
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I. INTRODUCTION 



Neutrino interactions in the low and middle energy region are attracting considerable 
attention because they investigate new properties [1, 2] and may lead to new discoveries. 
For this program precise estimates in the low and intermediate energy regions are necessary. 
Over the past few years many articles have been published, sometimes with conflicting 
results. The difi^erences among the theoretical conclusions and some of the resolutions have 
been summarized in a recent review article [3]. We mention three issues relevant to our 
work. 

1. There were two values for the dominant axial matrix element for the nucleon-A tran- 
sition. The PCAC prediction [4, 5] with the value C^(0) = 1.2 and a smaller value 
of 0.87 from a fit of the ANL data [6] with a large nonresonant background. After 
several studies which modified the background, the results of simultaneous fits for the 
ANL and BNL data give a larger value C^(0) = 1.10 ± 0.08, which is now accepted 
and prefered [7-9]. 

2. A sizable background under the delta resonance was introduced based on diagrams 
from chilal symmetry [6, 10]. Electroproduction data [11] gives a small background. 
Furthermore, apphcation of PCAC relates matrix elements of the axial current to 
pion-proton and pion-neutron data where a nonresonant background is very small. 

3. For experiments on nuclear targets a generator for the transport model must be intro- 
duced. The GiBUU code has difficulties in reproducing the MiniBooNE data [12-15]. 
Their predictions lie far below the experimental data. 

What we need now are reliable estimates and calculations. For this reason we went back 
to our earlier calculations [9] and selected the low region as a rehable benchmark for 
estimating cross sections. In this region the axial contribution is related by PCAC to cross 
sections and the vector contribution by CVC to data again. The introduction of experimental 
data includes nonresonant contributions. 

As another example of a discrepancy and its resolution, we mention the peak and the 
turn over in the differential cross section as — > 0. We found that in addition to the 
vanishing of form factors at the boundary of the kinematic region, the axial contribution 
has a peak at ^ 0.03 GeV^, refiecting the resonance peak in the integration over W. An 
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earlier result for the differential cross section was presented and is slightly improved in our 
figure 1, where the cross section is shown as the sum of various contributions. In figure 1 the 
prominence of the peak for the axial contribution is evident. A second property in this figure 
is the fact, that the vector-squared and the interference contributions are close to each other 
and add up for u^p — )■ /i~A++. For the antineutrino reaction V^n — )■ /i+A~ they cancel each 
other leaving the axial current contribution as the dominant one, with the smaller terms also 
included in our calculations. This introduces a new test for PCAC and, if confirmed, will 
give a reliable cross section for antineutrinos. The general trend is qualitatively confirmed 
by experimental results where the antineutrino reactions are always smaller [16-18]. 

In this article wc will extend the method to other reactions and will compare them with 
experimental data when available. We present the differential cross sections and we 
shall compare them with MiniBooNE results [19, 20] including nuclear corrections which 
are explained below. Since the method, including nuclear corrections, is transparent we 
can identify each contribution explicitly. This gives, in most cases, consistent results and 
whenever we encounter differences we discuss the modifications that may be necessary. The 
method of concentrating in the low region has been adopted by another group [21] which 
proposes to extend it to higher resonances by including a model with coupled channels for 
the production. 

Two features are special in neutrino reactions: 

1. the possibility to use various types of neutrino beams u^, u^, i/^, and 

2. the necessity to have various nuclear targets. 

Both bring peculiarities which we will address. The various neutrino beams will be present in 
long-baseline experiments because starting with a beam a component will be generated 
through oscillations. In the kinematic region of this article the variable is of the same 
magnitude as and m^. We use formulas with exact masses and include all helicitics for 
the lepton vertex. A difference is expected between i/^ and reactions which we calculate 
and present in section II. 

For nuclear corrections we use a model inspired by Fermi [22, 23], which is based on an 
analytic solution of the transport equation with the results summarized in a simple way, 
so that the various steps are evident. This allows one to identify the origin of the effects. 
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The direct inputs for the initial cross sections combined with the rescattering model brings 
a rather good agreement with the results of the MiniBooNE experiment [19, 20]. 

An earlier comparison [9] with the Argonne and Brookhaven experiments, in this limited 
kinematic region, was consistent. The early experiments have limited statistics for the z/^n — )■ 

A"*" channels. For this reason the determination of the nonresonant amplitude with / = 1/2 
is at best very approximate. Electroproduction experiments have a small J^l"^ amplitude [11]. 
Furthermore the reactions Tr^p — > Tr^p and the charge exchange 7r~p — > 7r°n are dominated, 
in this energy range, by the delta resonance, allowing only a small background [24, 25]. 

All results support the view that our approach is a good first order approximation. The 
fact that the nuclear corrections bring the theoretical curves closer to the data is well un- 
derstood and gratifying. 

In section II of this article we describe the method and present results for charged and 
neutral current reactions. In section III we apply the method to the MiniBooNE results and 
find good agreement. In the last section we present our conclusions. 

II. DESCRIPTION OF THE METHOD 

We consider the reactions listed in table 1. The reaction in the first row for A"^"^ pro- 
duction gives the larger cross section and has been studied most extensively. The Clebsch- 
Gordan coefficients (CGC) in the fourth column refer to the J^l"^ amplitude which is dom- 
inant and should give the prominent features of the reactions. In this article, however, we 
calculate each amplitude separately by using the following method. For the axial current 
alone we use the cross section 

(y{j:'^'p — >■ 'K^'p). (1) 

The matrix elements Loo, -^/zo ^-nd Lu were introduced in references [26, 27] and include 
the mass of the charged Icpton in the final state. The cross section (T(7r^p — )■ vr+p) is the 
production of hadrons at the center of mass energy W for the reaction i/^p ij,~p7r~^ and 
is taken from [24]. In the present calculation the Tr+p Tr+p cross sections are evaluated 
at W instead of u and this modifies slightly the curves in figure 1. A second improvement 
is the appearance of the additional terms proportional to L;o and Lu. For each reaction we 
use the appropriate experimental data. For instance, for the reaction i/^n — )■ /x~p7r° we use 
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FIG. 1: Anatomy of the various contributions to the cross section at E'i^ = 1 GeV. 

the reaction a{'K^n — >■ 7r°p) = a{'K~p — >■ 7r°n) from [25] as we shall discuss below. The fact 
that we use data means that the nonresonant background is already included. 

For the vector and interfererence terms we use the formalism of reference [4, 5], where the 
vector form factors where determined for the 1 = 3/2 amplitudes. Precise electroproducrion 
data established a small / = y2 nonresonant amplitude which we also include later by 
augmenting by 5% [10]. For the vector and interference terms we use the formalism of 
reference [4] with the vector form factors from [9]. For the vector-axial interference C^C^ 
we use the form factor C^{Q'^) extracted through PCAC, Cf — —jC^ and the vector form 
factors just described. The second column in table 1 gives the sign of the interference term 
which is constructive for neutrinos and destructive for antineutrinos. With these inputs we 
calculated the various contributions to the reaction z/^p — )■ fi~pn^ shown in figure 1. The 
curve denoted as rest indicates smaller contributions from additional form factors beyond 
CY and C^. 

We repeated the computation for the reaction i/^n — > jjTpi^^ using data from [25] . We shall 
use both cross sections in the comparison with the MiniBooNE data [19, 20]. The results are 

shown in figures 2(a) and 2(c). For antineutrino reactions the sign of the W3((5^, i^) changes 
as indicated in table I. The antineutrino differential cross sections are shown in figures 2(b) 
and 2(d). 

For an electron-neutrino beam we must introduce the very small mass of the electron. 
The result for electron- and muon-neutrino cross sections are shown in figure 3 for various 
incident energies. For comparison we include in the same figures the cross sections for 
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TABLE I: Input quantities and isospin factors for various reactions. 
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FIG. 2: Charged current differential cross sections for Ei, = 1 GeV. 
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FIG. 3: Comparison of muon and electron neutrino or antineutrino cross section for various energies 
-Ei^.F = 1, 2, 5, 10 GeV (bottom to top). The solid line is for electron neutrinos or antineutrinos and 
the dashed one for muon neutrinos and antineutrinos, respectively. 

beams. The muon mass turns the cross sections to zero as — > 0. The mass effect for the 
fgS is invisible. 

For neutral current reactions there are more changes. The effective interaction is 

^efr = (1 - 75) ^ [xVl + yAl + tVJ] (2) 

with and A^^ the isovector and V° the isoscalar hadronic currents. The parameters in the 
hadronic current are given in terms of the weak angle 9w 

2 

X = 1 — 2 sin^ y — —1 and 7=— -sin^^vK (3) 

with sin^ 9w ~ 0.25. The value ofy — —1 gives a constructive W3 interference term (because 
of the structure of the lepton current (1 — 75) u), making the neutrino reaction larger 
than the antineutrino. Beyond these parameters there is an overall normalization factor in 
the amplitudes. In the charged current interaction appears the current 

^i + aj = y2(:i^i±|5)=^^;r (4) 

and for the neutral current A^^. The CGCs are vahd for the triplet {A'j^ , A^^, A~) . Since 
we have taken the amplitude for the reaction 1/^^ — > /i~A++ as the standard amplitude, we 
divided the neutral current CGCs in the fourth column of table I by \/2. 
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FIG. 4: Neutral current reactions at = 1 GeV. 
An additional property of neutral current reactions is 

a{pp i/A+) = cr(i/n i/A°) 



(5) 
(6) 



which follows from isospin symmetry which for neutral current interactions is broken by the 
small term ^V^. The calculated differential cross sections for neutral currents are presented 
in figure 4. The equalities in equations (5) and (6) were also confirmed by the results of an 
analytical calculation in table 2 of reference [6] . The zero mass of the neutrino produces the 
finite values for the cross sections at = 0. This is also the exact point determined by 
PCAC where the neutrino and antineutrino cross sections are equal. 



III. COMPARISON WITH RECENT EXPERIMENTS 

In the resonance region pions are produced by charged and neutral currents. For our work 
the neutral current reactions are not relevant because the experiments do not determine the 
momentum transfer squared, since the momentum of the final neutrino is not observable. 

They measure the momentum spectrum of the produced pion, which we cannot calculate 
since it involves the integral over small and large values of Q^. For this reason we calculate 
the charged current differential cross section for tt"*" and 7r° production and compare them 
with data. 
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In the previous section we described the initial interactions on free protons and neutrons. 
In the nucleus the production and development of the resonance and of the decay products is 
influenced by the medium. The effects of the medium are included in a transport matrix [23, 
28] whose nature is determined by an absorption term and the interactions of the pions. For 
this purpose we introduce an inverse interaction length k and a charge exchange matrix. 

The inverse interaction length is assumed to be the same for all pions, given by 



with p(0) the nuclear density at the center of the nucleus and h±{W) describes the forward- 
and backward-hemisphere projections of the Pauli blocking factor and Uabs is the absorption 
cross section given in equation (27) of [23]; it is a phenomenological function which contains 
several effects: the excitation of the nucleus, the propagation of the delta resonance before 
decay etc. Some authors compute the interaction of the delta with the medium as a self- 
energy correction [29] which shifts the mass and width of the resonance. Introducing a 
Breit-Wigner form for the resonance (with a corrected width F -|- 5T) and taking its Fourier 
transform it becomes evident that a shift in the width corresponds to a shift in the mean 
free path 



with z — vt the length of propagation and v the velocity. This shift in the width of a 
resonance in momentum space corresponds to a shift in the inverse interaction length in 
configuration space. Thus using aabs accounts for several effects. 

During the propagation, scatterings take place as described by the pion-nucleon cross 
sections at the center of mass energy of the delta resonance. The final result M{(iC^^) has 
a simple form with an effective absorption factor A{Q^) and a charge exchange matrix. For 
the carbon target [28] 



(7) 



(w) [h+{w) + h4w)] 



(8) 





0.83 0.14 0.04 



M{eC'^) = A(Q2) 0.14 0.73 0.14 



(10) 



V' 



0.04 0.14 0.83 
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with 



A{Q^ = 0.05 GeV^) = 0.71 (11) 
A(Q^ = 0.20 GeV^) = 0.79 (12) 
A{Q^ = 0.40 GeV^) = 0.81. (13) 

The pion multiple-scattering is solved analytically as a stochastic problem [23] . The results 
of calculations confirm an old suggestion that charge- exchange corrections are substantial. 
They also follow a general principle: In a lepton-nucleus interaction the pions which have 
the same charge as the current are reduced. For pions with different charge than that of the 
exchanged current, the cross section is enhanced. 

We apply this formalism to the MiniBooNE results. The target in the experiment is the 
molecule CH2, which we consider as the incoherent sum of C^^ and two protons. This is 
justified since the two structures are relatively apart in the molecule. The final yields of 



and TT are indicated by S-* + and E'^o and are obtained from 

= A{Q') {0.83 + E:^] + O.UE^ojG + 2S^^ (14) 

E^o = A{Q') {0.73E;fo + 0.14 [E^^ + E;;+] } 6. (15) 

The cross sections within the brackets are defined as 

- ^ ^ ^^-P-^) (16) 

^ W ^ '""''''^^ ^^^^ 

= % {ly.n ^ t,-pn') (18) 



which we call the initial or primitive cross sections. The bars over the cross sections indi- 
cate averaging over the flux for 0.50 GeV < < 2.00 GeV and are the cross sections we 
computed for the spectrum of the MiniBooNE experiment. 

The final yields for tt"*" and 7r° are shown in figure 5 where we plotted also the data 
from MiniBooNE. As the vector contribution from does not include nonresonant back- 
ground [5] we augmented this form factor by 5% as a conservative estimate in this calculation. 
For the effective absorption factor A((5^) we use the values from equations (11)-(13). The 
agreement of the vr^-yield with the data is very good. For the 7r°-yield our curves are a little 
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(a)7r"'" yield (b)7r'^ yield 

FIG. 5: Charged and neutral pion production for the MiniBooNE spectrum and comparison with 
experimental data. The dotted line is the cross section without nuclear corrections, the dashed line 
with constant absorption A = 0.9 and the solid line with A{Q'^) interpolated from equations (10) 
to (12). 

lower. However, in view of the larger experimental errors and the simplicity of the model 
the agreement is satisfactory. 

On the experimental side, the data, especially the vr'^-yield, still have large error bars 
and on the theoretical side it may be necessary to add a background to the z/^n iJ,~A~^ 
production. Looking into the composition of the yields, one finds the tt"*" yield comes pri- 
marily from the direct production of tt'^s with a small feeding from the tt^s. The 7r° yield, 
on the other hand, receives a substantial feeding from the primitive tt"^. The contributions 
from the two terms, direct and feeding, are almost equal in this channel, as the reader can 
easily verify using the monoenergetic initial cross sections in figure 2 and the equations in 
this section. However, as we mentioned above, for the curves in figure 5 we averaged over 
the neutrino spectrum. The contribution from charge exchange (feeding from the Tr"*" yield) 
compensates the reduction from A{Q'^) and produces the solid curve in fligure 5. 

IV. CONCLUSIONS 

1. We presented a calculation for the production of pions by neutrinos in the small 

region and for low energies (0.50 - 2.00 GeV), based on PCAC, CVC and experimental 
data. We selected this region because the method is based on tested physical principles. 
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The results agree with the Argonne and Brookhaven data for the 7r+ channel, as has 
been shown earlier [9]. The results are consistent with a value for C^{Q'^) close to 1.1 
and a small background. A small background is also observed in electro-production [11] 
by comparing the 7r° and channels. A small background was also obtained in 
calculations at the static limit [30, 31]. The same Argonne experiment requires a large 
/ — 1/2 component, however one must keep in mind, that the statistics for the reaction 
f^n — >■ //~A+ are limited. 

2. We applied our method to the MiniBooNE experiment, which observed two different 
final states (Itt"*" and Itt^), and found good agreement, which further supports our 
estimates of the initial interactions and in addition for the transport model for the 
carbon nucleus [28]. For these results we used flux averaged cross sections. The tt^ 
cross section is affected by absorption and Fermi blocking with a very small correction 
from charge exchange. In contrast the 7r° channel receives a 45% contribution from 
charge exchange as the large ir'^p channel feeds into it. This is very different from the 
results of GiBUU [15] where they reported (see end of section III. B in [14]) that the 
"FSI have a relatively small effect" on the 7r° yield by reducing it, instead of increasing 
it. 

3. One of the uncertainties in neutrino reactions is the presence of many form factors, 
whose dependence is, in principle, unknown and determined from data. One can 
use for the leading form factor, especially for C^{Q'^), the functional form we found 
and match it to the extrapolation to larger vales of (see figure 3 in [9]). Thus we 
consider the results at < 0.20 GeV^ as a starting point for the parameterizations 
to larger values. 

4. The charge exchange matrix in equation (10) and the primary interactions obey isospin 
symmetry. Thus for isoscalar targets, like carbon, there are isospin relations between 
various reactions which can be tested in the new experiments (Minerva). In an internal 
report [32] we elaborated on some of them. They should be tested in the experiments. 
Finally, the input parameters for the initial reactions [9] and the nuclear corrections [28] 
were published earlier and we did not need to modify them. 
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